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Abstract 

f~^ ■ We investigate the behaviour on approach to the initial singularity in higher-order extensions of gen- 

O^ ' eral relativity by finding exact cosmological solutions for a wide class of models in which the Lagrangian 

is allowed to depend nonlinearly upon the three possible linear and quadratic scalars built from the Rie- 

mann tensor; R, RabR"'' and Rated R"^'"^'' ■ We present new anisotropic vacuum solutions analagous to the 

qH| Kasner solutions of general relativity and extend previous results to a much wider range of fourth order 

I . theories of gravity. We discuss the implications of these results for the behaviour of the more general 

w rv' anisotropic Bianchi type VIII and IX cosmologies as the initial singularity is approached . Furthermore, 

we also consider the existence conditions for some other simple anisotropic Bianchi I vacuum solutions 

in which the expansion in each direction is of exponential, rather than power-law behaviour and their 

relevance for cosmic "no-hair" theorems. 

> 

a^. 

}Q ■ 1 Introduction 

f^*- ' The literature contains numerous studies on generalisations of the familiar Einstein-Hilbert action of 

(«_^ I general relativity (GR) to more complicated functions of the curvature and higher-order invariants [HID 

>--^ ■ 13]. Motivation for these studies comes from several sources, including astronomical phenomena which are 

currently inadequately explained by the standard model of general relativity, such as providing a natural 
source of inflation in the early universe 2], or dark energy and the late-time acceleration of the universe's 
expansion [S] [SI [3111111 HO] > a^nd also attempts to include quantum behaviour in the gravitational theory 
[11) . A review of one of the most common extensions to general relativity, the so-called f{R) models, in 
which the Lagrangian is allowed to be a general function of the scalar curvature, may be found in [12) . 

It is of particular interest to discover the behaviour of these higher-order theories at high curvatures 
and it is in this limit when we might expect the influence of quantum corrections to become important. 
Therefore, where initial singularities are expected to involve infinities in one or more of the curvature 
invariants of the space-time, we expect that the addition of higher-order terms to the Lagrangian might 
produce a new dominant behaviour to such singularities. The (past) stability properties of special initial 
isotropic cosmological singularities were investigated in [IS] for higher-order theories where the dominant 
term in the Lagrangian took the form {RabR"''')"' ■ 

It is well known that contributions to the Lagrangian from terms dependent on the scalar curvature 
only are conformally equivalent to the presence of a minimally coupled scalar field in general relativity [TJ , 
however, in theories in which the gravitational Lagrangian contains higher-order curvature invariants, a 
much richer diversity of anisotropic behaviour is possible. For example, speciflc counterexamples were 
found in [TS] which demonstrate that the "cosmic no-hair theorem" of general relativity may be violated 
in higher-order theories. Furthermore, anisotropies diverge faster than isotropics at high curvatures and 
will tend to dominate the cosmological behaviour at early times. Thus, whilst the majority of previous 
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studies of these modified theories of gravity have focussed on the behaviour of isotropic cosmologies, it 
is the role of anisotropy on approach to the initial singularity which we wish to investigate here. 

Some of the most important anisotropic cosmological solutions in general relativity are the vacuum 
Kasner solutions of Bianchi type I 16 . Since they are characterised by just a single free constant, they 
are geometrically special, but nevertheless they provide us with a very useful insight into the dynamics of 
anisotropies, since they give a good description of the evolution of more general anisotropic cosmological 
models over finite time intervals. The chaotic oscillatory behaviour of the spacetime on approach to 
the initial singularity exhibited by the Bianchi type VIII and type IX ( "Mixmaster" ) cosmologies can 
be approximated by a sequence of different Kasner epochs [171 1181 1191 1201 1211 122| . Provided that at 
least one of the Kasner exponents is negative, inhomogeneities and perturbations from the Bianchi I 
anisotropies will grow as the singularity is approached and force the solution to switch from one set of 
Kasner exponents to another. If the solution must always have at least one negative Kasner exponent, 
as is the situation in general relativity, then these oscillations will continue infinitely as the singularity 
is approached. However, in some higher-order theories of gravity [2311241 E5] it may be possible for all of 
the Kasner indices to take positive values, whence after a sufficient (finite) number or permutations, the 
indices will reach such a configuration and the oscillatory behaviour will cease. Since all spatial directions 
will be contracting, the initial singularity will be reached monotonically. 

Previously, Kasner-type cosmological models in quadratic gravity and in Lovelock theories of grav- 
ity in higher dimensions were investigated by Deruelle [23] . Clifton and Barrow [241 125| discovered the 
conditions for the existence of Kasner-like solutions and the exact forms of these solutions for the par- 
ticular cases where the Lagrangian is an arbitrary power of one of the curvature invariants, R, RabR°'^ 
or RabcdR^'""^- One might expect that the dynamics and the asymptotic behaviour of any solution in 
a more general higher-order theory would be controlled by the highest powers of the curvature in the 
past, and the lowest powers of the curvature in the future. However, this assumption is not necessarily 
accurate [26] and we wish to extend the investigation of [241 125| to include more general Lagrangians. 

Thus, in this work, we wish to investigate, within this class of higher-order theories of gravity, the 
constraints on the Lagrangian for the existence of some simple anisotropic but homogeneous solutions 
of Bianchi Type I and to find all such solutions. In particular, our main focus will be to consider the 
possibility of anisotropic Kasner-like solutions in vacuum, and in the presence of a cosmological constant. 
We will also consider the properties of these solutions with respect to their relation to the behaviour of 
the more general Bianchi type VIII and IX cosmologies. In addition, we will also discover all solutions 
in these higher-order theories which are expanding anisotropically but exponentially in each of the three 
spatial directions. 



2 Field Equations 



In this paper we will consider theories of gravitation in which the field equations are derived from 
an arbitrary analytic function of the three possible linear and quadratic contractions of the Riemann 
curvature tensor; R,RabR°''' and RatcdR'^'"^'^ ■ The relevant action is given by 

-f{X,Y,Z) + L„ 

A 

where f{X, Y, Z) is an arbitrary function of X, Y and Z which are defined X = 7?, Y = R'^^Rab and 
Z = R'^^'^'^Rabcd- Spacetime indices run from to 3 and are denoted by roman letters, whilst Greek 
letters are used to denote purely spatial indices. 

Unlike the situation in general relativity, which may be recovered by choosing f = R, the Palatini 
and metric formalisms are not equivalent for a general choice of /. In what follows, we shall restrict 
attention to the metric formalism. The field equations obtained by varying the action with respect to 
the metric are [25[ : 

Pb = fn" , (1) 



where 



Pab -L /■ ah , r ■n>o,b , r» r Ticfa i~,b) , r» r -medcia i~,b) , / ab cd ac bd\ c 

= -0/5 +/X-R +2fYR' R' c + 2fzR 'R'cde + {g g - 9 g )l 
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+ n{fYR'"') + g"-\fYR''');cd - 2(fYR''-''); e' ^ 4(/zi?'<"''>") ;cd : 



X;cd 



where fx = -q^, fv = -gy ^nd fz = g^- Whilst the focus of this work is on solutions in vacuum, 
with Tl^ ~ 0, the study can be extended in a simple way to allow the possibility of a non-zero cosmo- 
logical constant, A, by including any such term in the gravitational part of the Lagrangian, f{X, Y, Z). 
Furthermore, we will also consider the effects of including a perfect fluid for f{R) theories of gravity. 

3 The Kasner Model 

In this paper, we consider homogeneous but anisotropic Bianchi I models, and in particular our main 
focus will be on those spacetimes described by the line element 

ds^ = -dt" + t^^^dx' + t'^'^dy' + t'^'-'dz^, (2) 

where the Kasner exponents pa are constants and assumed to be real in order for the metric to be of 
physical significance. We define the useful quantities 

H = Pl + P2 + P3 , 

J = Pi +P2 +P3 and 

T^ — 3 , 3 , 3 

K = Pl +P2 +P3 

The line element ([2} also describes the limiting case in which all three Kasner exponents are equal. 
Pi = P2 = P3, and corresponds to isotropic, spatially flat Friedmann- Robertson- Walker solutions for 
which the scalefactor is a power-law in time. These solutions will be included for completeness. It is 
useful to note that for real- valued choices of pa , J > ■%- , with equality if and only if the solution is 
isotropic. The three relevant curvature scalars with which we will be working take the values 

Y^R.,R-' = AH-lf + iJ^Hf ^ ^^^ 

7 _ „ pabcd (3J-H^)^ + 12J + 8J^(H-3) ... 

Z = liabcdil — r-j • (O) 

In general relativity, the Kasner exponents for the vacuum solution must satisfy H = J — 1 |16| . 
One solution of this is the Milne model [27] , which without loss of generality can be described using the 
choice of axes such that pi = 1, P2 = Ps = 0. However, this solution is related to Minkowski space by 
a coordinate transformation; if we introduce new coordinates r — icoshx,x ~ isinha;, then the usual 
form of the flat Minkowski metric is explicitly recovered [25] • The Riemann tensor Rated vanishes for 
Minkowski space, and therefore this is a vacuum solution in any higher-order theory of gravity of the 
form f{X, Y, Z) for which /(0, 0, 0) = 0. 

For all other general relativistic Kasner solutions, one Kasner exponent must be negative, whilst the 
other two are positive. Thus, although the spacetime volume is expanding to the future, one of the 
spatial directions is contracting. Moreover, as the initial singularity is approached, inhomogeneities and 
deviations from the Bianchi type I anisotropics in the more general Bianchi type IX solution will grow 
and cause the Kasner exponents to be permuted to different values, leading to an infinite series of chaotic 
BKL osciUations between different Kasner epochs [171 [TSl [T9l [20l [211 [22] . 

However, this is not true in general for the higher-order theories considered in this paper. In [241I25| . 
Clifton and Barrow found some exact solutions, a subset of which permits all the Kasner exponents to 
take positive values. We will see that this such solutions exist in a much wider class of higher-order 
theories. In this scenario, it is possible that after a finite number of transitions between different Kasner 
epochs, the solution will reach a state in which all of the Kasner exponents are positive. Once this occurs, 
the perturbations from the Bianchi I model would not grow and the chaotic oscillations will cease. Thus, 
the solution will remain in this epoch and the initial singularity is then approached monotonically. 

By considering those classes of theories in which only one curvature scalar - $, say - contributes a 
time scale to the Lagrangian, one can solve for the time coordinate t in terms of the scalar $. In this 
way, time derivatives may be eliminated and the field equations may be re-written in terms of $, the 
Lagrangian /($), and its derivatives with respect to "I>. The resultant differential equation(s) can then 
be solved to find all possible forms of /. This technique, which was also used by Dunsby et al. in 



f{R) theories [28], allows us to find all possible exact Kasner-like vacuum solutions within this general 
class of Lagrangians, with one exception, which may be dealt with separately. This exception is when 
quantity which usually determines the time scale becomes independent of time for some special choice 
of the parameters pc,. For example, in general for the metric ([2}, the scalar curvature, R, is proportional 
to t~'^, and so in the context of f{R) models one can substitute for the time coordinate t using R, but 
it is necessary to consider separately those possible solutions with R = 0. 

As a consequence, we will find that the subset of metrics for which the curvature scalars are indepen- 
dent of time plays an important role and it is useful to discuss those metrics briefly now. Note from the 
expressions ((3]), (|4]and ((5]), the scalars R, RabR"'^ and RabcdR"'^'^'^ are only constants if and when they are 
zero. A special case within this class of metrics is flat Minkowski space, which may be respresented by the 
metric ^ with pa = 0, that is to say gab = flab, for which the Riemann tensor, Rabcd, is identically zero, 
and thus X = Y = Z = 0. This is a solution in all theories for which /(O, 0, 0) = and the discussion in 
subsequent sections will be concentrated on metrics describing curved spacetimes. 

By expressing the Kretschmann scalar, Z = Rabcd R"^'"^'', as 

Z ^ — [pi (pi - 1) +p2iP2 - 1) +P3 (P3 - 1) +Pl P2 +Pl pa +P2 P3 ) , 

it can be seen that there are two possible real solutions to Z — 0, given by Minkowski space, pa — 0, 
and pi = 1,P2 ~ P3 = (plus permutations). However, as we have seen, the latter solution is the Milne 
model, and is related to the former by the coordinate transformation r — i cosh a;, x ~ isinha;. In order 
to satisfy Y = RabcdR"'^'^'^ ~ 0, real solutions must have either H = J — or H = J— 1. The only 
possibility in the former case is Minkowski space. The latter are the well-known Kasner solutions of 
general relativity [IS]. Thus, Z — implies that Y = 0. 

Solutions for which the scalar curvature, X = R, is zero require J = 2H — H^ and thus Y = implies 
that X = 0. From this equation for J, one can see that real solutions must satisfy the constraint J < 1. 
Furthermore, from our definitions, J > -^ for all real- valued choices of the constants Pa, and so it is 
also necessary that < H < 3/2. According to these constraints, we see that the spacetime volume 
of any solution of this form must expand no faster than i^", and that the expansion rate in any one 
direction can be no faster than t. Finally, we note that solutions with X = and H < 1 cannot have 
all three Kasner exponents positive and, except for Minkowski space and the Milne model, at least one 
exponent must be negative. This is the situation in general relativity and, except for the two cases above, 
implies that although the spacetime volume is expanding, space must be contracting in one direction. In 
contrast, if X = but 1 < H < 3/2, then it is possible for all the Kasner exponents to be distinct and 
positive. Thus, as we have discussed, these solutions can avoid the infinite series of chaotic oscillations 
between different Kasner epochs on approach to the initial cosmological singularity seen in the BKL 
picture. The only solution with X — and H = 3/2 is the isotropic solution with pi = P2 = Pa = 1/2. 

In what follows, we shall consider in turn several commonly-studied classes of higher-order theories 
of gravity, obtaining the forms of the Lagrangian within these classes for which solutions of the form ([2|> 
exist, and the conditions which the Kasner exponents are subject to in each case. 

3.1 f = f{R) 

If the Lagrangian depends on the scalar curvature only, / = f{R), then it is well-known to be conformally 
equivalent to general relativity with a minimally coupled scalar field 14 . Thus, one does not expect 
a large range of anisotropic behaviour to be possible in such models; indeed for the case of quadratic 
corrections to the Einstein-Hilbert Lagrangian it is precisely the presence of the Ricci term, RabR"^^, 
which permits the existence of the anisotropic Bianchi I solutions found by Barrow and Hervik [15] . 
The vacuum field equations are Pj" — 0, where, for the Kasner-like metric given by ^, we have 

po t^f + 2(H ^ J)fR + 2HtfR 

Po - ^2 ' (6) 

pc _ ~3t^f + 2H{H - l)/fl - 4mfR - et^fR 
° - 2i2 ' •■ '' 

PH-P: = ^{p^-p.){iH^l)fR+tfR) and (8) 

~f2 



« . i^igm^ (9) 



where ^, v in equation ((8)l for the anistropic stress are not indices to be summed over, but instead are 
used to label the Kasner exponents, taking values from 1 to 3. Otherwise summation convention is used 
as normal. Recall that we may allow / to contain a cosmological constant term, /ij is used to denote ^ 
and overdots represent derivatives with respect to the coordinate time, f. 

Given the form of equation ([U , in what follows we shall consider separately the situations where the 
metric is isotropic and where it is anisotropic, both for vacuum and with the inclusion of a perfect fluid. 



3.1.1 Isotropic power-law vacuum solutions 

If all the Kasner exponents are equal, the solution is isotropic, with pi 
and there is only one independent field equation: 

3tV - 2H{H - 3)fR + 6HtfR = . 



P2 = P3 



^ - 3l^ 



(10) 



For a particular theory, that is to say a particular choice of f{R), this equation might appear to provide 
a constraint on H, the one free constant remaining. However, in general, it is algebraic in H and t 
and therefore one will not always be able to find constant solutions of this equation for H; in fact such 
solutions are rare. 

For models with Lagrangians of the form / = f{R), it is possible to summarise the full set of possible 
isotropic power-law vacuum solutions, and the conditions on the function / for these to be valid, using 
the classification in table [T] 



Class 


Solution 


Validity 





H^O 


Minkowski space is a vacuum solution in any model with /(O) — 0. 


I 


J/ = 3/2 


This solution, analagous to the radiation-dominated Friedmann uni- 
verse of General Relativity, is a solution of the vacuum field equa- 
tions for any model satisfying /(O) — and /i{(0) = 0. 


II 


tr 3(2n-l)(n-l) 
-" 2-n 


This class of solutions is possible if and only if the Lagrangian takes 
the form / = Q„_R" + a,„7?™, where n ^ {0, 5/4, 2}, m = ^j^ and 
a„,am are constants, with q:„ 7^ 0. These solutions are expanding 
to the future if n < 1/2 or 1 < n < 2. 



Table 1: Isotropic power-law vacuum solutions in f{R) gravity. 



The class II solutions were found previously in [24l in the context of models where the Lagrangian is 
a power of the scalar curvature, / — R". However, it is important to note that for all such models with 
n > 1, the solution of class I, which has zero scalar curvature, also solves the vacuum field equations, a 
fact which was not stressed in [24| . Furthermore, for any model in which the Lagrangian may be written 
as a power series in R, a non-zero linear term (i.e. an Einstein-Hilbert term) in the series precludes the 
possibility of isotropic power-law vacuum solutions other than Minkowski space. 



3.1.2 Isotropic power-law solutions with a perfect fluid 

Let us now consider a universe filled with a comoving perfect fiuid, with equation of state p = wp, and 
energy density evolving as p{t) = poi"^ , where po is a constant. The possible isotropic power-law 

solutions of the field equations for a Lagrangian of the form / = f{R) are similar to those found in the 
vacuum case and are summarised in table [S] 

The solution of class I is analagous to the radiation-dominated Friedmann universe of General Rel- 
ativity, and also to the class I vacuum solution when an extra Einstein-Hilbert term is added to the 
Lagrangian. 

It is apparent that the vacuum solutions found in the previous section correspond to fiuid-fiUed 
solutions with an appropriate extra curvature term in the Lagrangian. 



Class 



Solution 



Equation of state 



Validity 



I 
II 

III 



H = ?,/2 



l-\-w 



This is a solution in any model satisfying 
/(O) = and /^(O) = constant / 0. 

This class of solution exists iff the La- 
grangian is a power of the scalar curva- 
ture, / — OnR", where n 7^ and a„ is 
a constant. 

This is a solution iff the Lagrangian takes 
the form / = a„R" + a^R"' + aj,RP , where 
n i {0, 5/4, 2}, m = ^i^^s. and a^ / 0, a„ 
and am are constants. For power-law La- 
grangians [am = ctn =0), this solution 
reduces to a subset of the class II solutions. 



Table 2: Isotropic power-law solutions in f{R) gravity with a comoving perfect fluid. 



H 



3(2n-l)(n-l) 
2-n 



W= 1/3 
W = -1 + 



2p(2 — n) 
3(2n-l)(n-l) 



3.1.3 Anisotropic vacuum solutions 

We have seen that if the Lagrangian is a function of the scalar curvature only, / — f(R), there is a strong 
constraint on its form for the existence of exact isotropic power-law solutions. We shall now turn our 
attention to the situation for anisotropic Kasner-like solutions within these models, the main focus of 
this paper. In this case, the vacuum field equations can be reduced to 



iH^l)fR+tfR = 0. 

/ = 0. 
RfR = . 



(11) 
(12) 

(13) 



If the arbitrary function / is non-trivial and algebraic in the scalar curvature, R, then the only 
possible solutions for R of f{R) = must be constants, and therefore any such solution satisfies R = 
and thus J = 2H — H^. Recall that from our definitions, solutions with zero scalar curvature for which 
the constants Pa take real values must satisfy the constraints < H < 3/2 and J < 1. According to 
these constraints, we see that the spacetime volume of any solution of this form must expand no faster 
than t''", and that the expansion rate in any one direction can be no faster than t. 

Whilst /(O) = is a necessary condition for a f{R) model to contain anisotropic Kasner-like solutions 
in vacuum, and this also implies that equation (|13p is satisfied at i? = 0, it is not sufficient, since /fl(0) 
is model-dependent. Thus it is not guaranteed that there will be solutions of equation (|lip with constant 
H in models where /ij(0) diverges; in this case, solutions will only exist if the divergent part of the 
Lagrangian is a power of the scalar curvature. If, however, /fl(0) is zero, then equation (|lip is satisfied 
trivially. If /ii(0) is a non-zero constant, as is the case in general relativity, defined by f{R) = R, then 
this gives the extra constraint H = 1, whence J = 1 and the only possible solutions are those of general 
relativity. 

To exemplify the situation, we could consider a Lagrangian which can be expanded as a power series 
about R = 0. A non-zero linear term - an effective Einstein-Hilbert term - in the series precludes 
the possibility of anisotropic Kasner-like vacuum solutions other than the general relativistic one, with 
H = J — 1, whilst a non-zero constant term - an effective cosmological constant - would preclude the 
possibility of these solutions altogether. 

A summary of the conditions that must be satisfied by the model in order for solutions of this kind 
to exist and classification of the possible solutions is given in table |3] 

The first class of solutions contains a subset for which all the Kasner indices may be positive; this 
is possible if 1 < H < 3/2. We have seen that the existence of such solutions means that the infinite 
series of chaotic oscillations on approach to the initial singularity in the more general Bianchi type IX 
Mixmaster universe can be avoided. 

The third class of these solutions are relevant for models in which the Lagrangian is an arbitrary power 
of the scalar curvature, / — aR", and were found in that context by Clifton and Barrow [24]. However, 



Class 



Solution constraints 



Validity 



I 
II 

III 



J ^2H-H^ 



H = J =1 



H = 2n-l, 

J = (2n-l)(3-2n) 



These solutions are possible if /(O) = and /i{(0) = 0. The 
exponents pa are real provided < H < 3/2. 

If the model satisfies /(O) — 0, but /fl(0) is a non-zero constant, 
then the extra constraint from equation (lllf) means that only 
this subset of the first class of solutions is possible. These are 
the solutions for the case of general relativity, f{R) — R. 

This subset of the first class of solutions is relevant if /(O) — 0, 
but ,ffl(0) diverges due to a term which is a power of the scalar 
curvature. Thus, the Lagrangian is required to be of the form 
/ = aR"+f{R), with < n< 1 and /(O) = /h(0) = 0, although 
/ need not be identically zero. These solutions are complex and 
of limited physical interest if n < 1/2. 



Table 3: Kasner-like vacuum solutions in f{R) gravity. 



for n > 1, we have found here that there are additional exact Kasner-like solutions corresponding to 
R = which do not necessitate that _ff = 2n — 1 and are subject only to J = 2H — H^. 

3.1.4 Anisotropic solutions with perfect fluid 

It is interesting to also include the possibility of a comoving perfect fluid with a barotropic equation of 
state, p — wp. The energy density of the fluid is given by 

-H(l+w) 



P = pot 



(14) 



If the fluid is comoving, the spatial part of the energy-momentum tensor is isotropic and so equation 
(|lip still holds. This can be solved to give /ij oc t^~^ and the remaining field equations simplify to 



p(l + w) 
■wp 



-RfR, 

-//2. 



(15) 
(16) 



The only solutions of these equations with constant non-zero energy density correspond to a cosmological 
constant, which was considered together with the vacuum case in the previous section. For all other fluid- 
filled solutions, p oc ^-^(1+™) and R oc t~^, and thus the field equations can only have a solution of this 
sort if the Lagrangian is a power of the scalar curvature, f{R) = -R", for some n £ R\{0}, and R ^ 0. 
For / = R^'^, the right hand sides of the above equations are identically equal, hence to allow a non-zero 
energy density, we further require n 7^ 1/2. 

For all n 7^ 1/2, these equations have the solution H = 2n— l,w = (2n — 1)~^. Since these solutions 
need 7? 7^ 0, it is also required that J 7^ (2n — 1)(3 — 2n), but otherwise it is unrestricted. 



3.1.5 Examples of solutions for specific choices of f{R) 

To summarise the results of the previous sections, in table |4] we consider some of the more commonly- 
studied choices for the Lagrangian, f{R), and state whether these models contain Minkowski space, 
isotropic power-law solutions, and exact anisotropic Kasner-like solutions, both in vacuum and in the 
presence of a comoving perfect fiuid, according to the classifications given in the preceeding sections. We 
include the choices of power-law Lagrangians [24] . and two exact f{R) models which have recently been 
proposed, by Starobinsky [5U], and by Hu and Sawicky ^D]. These models are of particular interest since 
they provide viable cosmologies and evade the known constraints on the form of f{R) from solar system 
tests. They are given by 



f{R) = fst.r{R) 
f{R) = fnu{R) 



R-\Ro[l- 
R + \Ro 



1 + {R/RoY 
{R/RoT 



and 



1 + a{R/RoY 



(17) 
(18) 
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respectively. The parameter n is taken to be positive for both the Starobinsky and the Hu-Sawicky 
models, in order to ensure their viability with observations. 



fiR) 


Minkowski 


Power-law 


Kasner 






in vacuum 


with fluid 


in vacuum 


with fluid 


R 


/ 


X 


Class I, II 


Class II 


/ 


R + A 


X 


X 


X 


X 


X 


R + aR^ 


/ 


X 


Class I,III 


Class II 


X 


R + a/R 


X 


X 


Class III 


X 


X 


R",n<Q 


X 


X 


Class II 


X 


/ 


R",Q<n<^ 


/ 


Class II 


Class II 


X 


/ 


R'/^ 


/ 


X 


X 


X 


X 


R",\<n<l 


/ 


Class II 


Class II 


Class III 


/ 


i?",n> 1 


/ 


Class I, II 


Class II 


Class I 


/ 


exp (R/Ro) 


X 


X 


X 


X 


X 


/star 


/ 


X 


Class I 


Class II 


X 


/Hu,n> 1 


/ 


X 


Class I 


Class II 


X 


/hu, 1 > n > 


/ 


X 


X 


X 


X 


/Hu,n< 


X 


X 


X 


X 


X 



Table 4: The existence of power-law and Kasner-likc solutions in various models of f{R) gravity. 



3.1.6 Remarks 

Friedmann- Robertson- Walker power-law solutions in f[R) gravity with a perfect fluid were investigated 
in [28]. There, it was claimed that the only possible form of Lagrangian which admits these solutions 
and has the correct general relativistic limit is a power-law, / = R" . In fact, there is another class 
of isotropic power-law solution (type I in our classification) corresponding to a radiation-filled universe, 
which is possible in any theory for which /(O) = and /fl(0) is constant, however comoving perfect fluids 
with other equations of state are not possible. Furthermore, we have seen here that it is only the 7?" 
Lagrangians which allow anisotropic Kasner-like solutions with a perfect fluid. We can see, therefore, 
that these i?" theories of gravity are special in admitting this sort of solution. Some Bianchi type I, III 
and Kantowski-Sachs solutions in f{R) gravity have also been investigated recently by Farasat Shamir 

m- 

Anisotropic singularities of Bianchi Type I were recently studied in the context of more general 
Lagrangians of the type f{R,(f),x), with x = —^g"^''da4>db4> [32], with particular focus on the anistropic 
instabilities related to the existence of the hypersurface 4^ — and solutions being able to cross this 
surface, leading to questions about the viability of these models. Here, it has been shown that exact 
anisotropic Kasner-like solutions in f{R) theories must have R — and although they can live on the 
hypersurface defined by g^ = 0, they cannot cross it. 



3.2 f = f{R-^Ra,) 

Let us now consider the case where the Lagrangian is a function of the Ricci invariant only. For isotropic 
cosmologies, the contributions to the field equations from the simplest such term, RabR'^'', are propor- 
tional to those from a term quadratic in the scalar curvature, R^. However, this is not true more generally, 
and the Ricci term allows much more diverse anisotropic behaviour [15] . 

For this class of theories, the relevant field equations for the metric ([2]) in vacuum are P^ = 0, where 






1 

2t4 
1 



-t'^f + 2{-H^ 



H' 



HJ + 2/ - J)fY + 2{H'' + J- 2HJ)tfY 



(19) 



(p^ - p.) 2{H - 3)(J - l)/y + (-4 + 2J + 3H 



H^)tfY + {l-H)t\fY] , (20) 



where as before /.i, u in equation (|20|l for the anistropic stress are not indices to be summed over, but are 
used here as labels, taking values from 1 to 3. Otherwise summation convention is used as normal. Recall 



that fy is used to denote ^ and overdots represent derivatives with respect to the time coordinate t. 
The Ricci term is given by 



Y : 



, J{H-lY + {J- 



Hf 



f4 



(21) 



It is useful to recall our earlier observation that for real- valued choices of the constants pa , the Ricci 
term, Y, takes non-negative values, and is zero if and only if// = J=lorii' = J = 0. 



3.2.1 Isotropic power-law solutions 

For / — f{Y), and the special case of an isotropic metric, pi = P2 = Ps = H/3, so that J = H^ /3, the 
equations reduce to: 



9iV + 4ff^((-6 + 3H + H'^)fY - 3{H - 2)i/y) = . 



(22) 



By an argument analagous to that used before for the f{R) theories, f{Y) can be zero for all times only if 
y is a constant, and therefore zero. Using (|2H) . it is clear that Minkowski space is the only real isotropic 
power-law solution for which Y — 0. If y is non-zero, then one can eliminate the time variable t and 
instead consider this equation as a differential equation in Y . It may then be integrated to find that, in 
order to admit solutions of this sort, the Lagrangian is required to be of the form of a power of the Ricci 
invariant, or possibly a sum of two such terms. 

We can now summarise the existence conditions for vacuum solutions of this type for theories with La- 

grangians of the form / = f{Y), using the definitions H±(m) 



and n±(m) 



+ 



) 2 + 10m-24m2=|:2-^3(-l + 10m-5m2-40m3+48m'l) 



3-9m + 12m^±-^3(-l + 10m-5m2_40m3+48m'l) 
2(l-m) 

This is found in table [S] 



Class 



Solution 



Validity 




I 

Ila 



lib 



H = Minkowski space is a solution in any model with /(O) — 0. 

H = 3/2 This is a solution if the Lagrangian is linear in Y, f{Y) — aY, 

where a is a non-zero constant. 

H = H+im) This solution is possible if and only if the Lagrangian takes the 
form / = ctmY'^ -f a„y", where m, 7^ 1, n = n+{m and a„, am 
are constants, with am 7^ 0. These solutions are expanding to 
the future if tti < 1 

H = H^ (m) This solution is possible if and only if the Lagrangian takes the 
form / — UmY™' + a^Y" , where m j^ 1, n — n-{m) and an,am 
are constants, with am 7^ 0. These solutions are expanding to 
the future if m < 1/4 or m > 1/2. 

Table 5: Isotropic power-law vacuum solutions in /(_/?„(,_/?"'') gravity. 



As with the situation for the Lagrangian / = i?^, we see that the radiation-dominated Friedmann 
universe of General Relativity is a vacuum solution of the higher-order quadratic theory. 

Whilst H+ is unbounded as n — > ±cx3, H- is bounded and tends to 2 in both these limits. For 
power-law Lagrangians in the Ricci term, both class Ila and lib solutions are valid; these were found in 
[25] and their stability on approach to the initial singularity under small perturbations of the metric was 
previously studied in [13) . 



3.2.2 Anisotropic vacuum solutions 

For anisotropic Kasner-like solutions, we again consider separately the cases Y — and y 7^ 0, so that in 
the latter scenario, we can replace the time variable i and treat the field equations p9j|20p as differential 
equations in the Ricci term, Y = RabR"^^- Similarly to the situation for isotropic power-law metrics in 
these theories, anisotropic solutions of the field equations with y 7^ can only exist if the Lagrangian 
for the theory is a power of the Ricci term, with / = aY". 



Class 


Solution constraints 




Validity 


I 


H = J 






There is a family of anisotropic solutions of this sort if 
f{Y) = aY^'^ . The Kasner exponents pa are real pro- 
vided that < i:f < 3. 


II 


H^J=1 






These solutions require that /(O) = 0, and also that 
/y(0) either converges to a constant or diverges slower 
than y-i/^ 


III 


H^ (l-2n)2, 
J = (l-2n)(l- 


- 6n 


+ 4n2) 


These are solutions for Lagrangians of the form f(Y) = 
/3Y" . In order for the Kasner exponents to be real, it is 
needed that (1 - 2n)(l - 6n + 4n^) > 0. 



Table 6: Kasncr-like vacuum solutions in f{RabR"^^) gravity. 



The conditions for the existence of anisotropic Kasner-like solutions within this class of theories and 
the constraints that must be satisfied by the Kasner indices may be summarised as in table |6l 

The first class of solutions have real Kasner exponents provided < -ff < 3. The volume of 
these solutions can thus expand no faster than t^, but the individual exponents must lie in the range 
I (l — \/3) < Pa < I (l + \/3) , and so the expansion may be faster than the speed of light in a particular 
direction. 

The second class of solutions are the general relativistic solutions, with H — J = 1; we have seen that 
these are the only anisotropic metrics of Kasner type for which the Ricci term is zero. It is interesting 
to compare these with the class II anisotropic vacuum solutions found in the previous section for f{R) 
models, which also have H = J = 1. In that context the solutions require that /ij(0) is constant, and 
thus / ~ i? as 7? — >■ 0, whilst here / ~ y^' ^ (or higher powers) as y — >■ 0. 

The third class are the solutions found in 25 , in which theories where the Lagrangian is a power 
of the Ricci term were previously investigated. The condition that the third class of solutions have real 
exponents is satisfied if either ni < n < n2 or 1/2 < n < n^, where ni,n2 and na are the roots of 
(1 — 6n + 4n^) = 0, chosen such that rii < n2 < n-j,. These roots have approximate numerical values 
n\ ~ —1.30,712 ~ 0.17 and nz ~ 1.13. These boundaries correspond to isotropic solutions and one can 
see their relevance for the existence of isotropic solutions in [13]. For n > 0, both H and J must be less 
than unity, and so for positive n the expansion of the solution cannot accelerate. However, H > 1 for 
these solutions if n < 0, and in particular for models with ni < n < ^ (l — \/3), H will be greater than 
3 and so the volume of the spacetime must increase faster than t'^ and undergo an accelerated expansion. 
Whilst this third class of solutions includes the special case of quadratic gravity defined by n = 1, it is 
easy to see that the conditions give the second class of solutions, H = J — 1, and there is no additional 
set of solutions in this case. 

3.2.3 Examples for specific f{R''^Rab) 

We are now able to summarise the possible vacuum Kasner-like solutions for various choices of Lagrangian 
of the form f(Y). In table[71 we consider several more common examples of this type of model and detail 
whether they allow isotropic power-law or Kasner-like solutions in vacuum, according to the classification 
systems we have used in the preceeding sections. In this table, use is also made of the definitions given 
in the previous section; recall that ni,n2 and na are the roots of (1 — 6n -j- An?) = 0, chosen such that 
n\ < n2 < n^, and have approximate numerical values ni ~ —1.30, n2 ~ 0.17 and 713 ~ 1.13. 
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fiY) 


Minkowski 


Power-law 


Anisotropic Kasner 


Y 


/ 


Class I 


Class II 


Y + A 


X 


X 


X 


Y",n<ni 


X 


Class II 


X 


Y",ni <n<0 


X 


Class II 


Class III 


Y",0 <n<n2 


/ 


Class II 


Class II,III 


F",n2 <n< i 


/ 


Class II 


Class II 


yi/2 


/ 


Class II 


Class I 


y",i <n<n3,rz/l 


/ 


Class II 


Class II, III 


Y",n3<n 


/ 


Class II 


Class II 


Y + aY^ 


/ 


X 


Class II 


exp{Y/Yo) 


X 


X 


X 


sin{Y/Yo) 


/ 


X 


Class II 



Table 7: The existence of power-law and Kasner-like vacuum solutions in various models of /{RabR"^^) 
theories of gravity. 



3.3 f = f{Z) 

For functions of the Kretschmann scalar, Z = RabcdR"^'"''' , only, the field equations for the metric © in 
vacuum reduce to Pj" = 0, where 



P" — 
P" = 

Z = 



-i/ + 2Zfz + ^{J-K) [{H - 3)fz + tfz 
2 ^P^-P"^ (2{H - p^- p^)t({H ^2)fz +tfz) + {Z- H){A- 



(23) 



f* 



2H + H^~ 3J)fz 



+(-8 + 8H- 3H^ + 3J)tfz + 2(1 

(3J - H^f + 12J + 8K{H - 3) 
3t* 



H)t'fz 



with 



(24) 



(25) 



Once again fi, v in equation (|24[) for the anisotropic stress are not indices to be summed, but are used 
to label the Kasner exponents, taking values from 1 to 3. Otherwise summation convention is used as 
normal. Recall that jz denotes ^ and overdots are used to represent derivatives with respect to the 
time coordinate t. 

Recall also that, by expressing Z explicitly in terms of the Kasner exponents, it can be seen that 
the only possible real solutions of Z = are Minkowski space and the Milne model, which may be 
defined by pi = l,p2 = Pa = and is related to the Minkowski solution by the coordinate transformation 
T = icosha;,x = isinha;. Thus, for all other Kasner-like metrics with Z 7^ 0, it is possible as before 
to substitute for the time variable t using Z and integrate to find the requirements on the Lagrangian 
/ = /(Z) in order to permit such metrics as solutions of the field equations. Using this method, it is 
found that in order for Kasner-like solutions to exist in vacuum, the Lagrangian must be a power of the 
scalar Z, or possibly a sum of two terms if they permit the same values of the Kasner exponents pa- 
Such power-law Lagrangians were first investigated in |25) . 

3.3.1 Anisotropic solutions 

For the previous two types of models that were considered, with / = /(i?) or / = f{R°'''Rab), there 
was only one independent equation for the anisotropic stress. In contrast, for models with Lagrangians 
that are functions of the Kretschmann scalar, / — f {R"''"^'' Rabcd) , there are in general two independent 
equations for the anisotropic stress. This appears to be due to the fact that the Kretschmann scalar 
cannot be expressed solely in terms of the parameters H and J, unlike the scalar curvature and Ricci 
scalar. For the cases of interest, where the Lagrangian is specified by / = aZ" , with n 7^ 0, these simplify 
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to 

= {pi~p2){l + H~4n){-4-3J + H{Q + H-8n) + 8n + 8p3{n-l)), (26) 

= {p2-pa){l + H-4n){-4-3J + H{6 + H-8n) + 8n + 8pi{n-l)). (27) 

Together, these imply that 

0= {Pi-P2){pi~pa){p2-Ps){l + H~4n){n-1) , (28) 

and so fully anisotropic solutions with no two Kasner indices equal can be possible only if either n = 1 
01 H = 4n~ 1. However, substituting H — 4n— 1 into the remaining independent field equation Pi = 
leads to the further requirement that Z = 0, for which we have seen that the only real solutions are 
Minkowski spacetime and the Milne model. The only other possibility for a fully anisotropic solution 
with Pa all real and distinct is in the theory with n = 1, i.e. / — aZ, a particular case of the commonly 
studied class of quadratic gravity theories, which will be discussed in more generality in the next section. 
Using the remaining independent field equation, such anisotropic solutions must have H — J — 1. Recall 
that these were also solutions for the other special cases of quadratic gravity considered in the previous 
sections; / = aR^ and / = aY. 

3.3.2 Locally rotationally symmetric solutions 

Unlike the situation for the previous two classes of model, with / — f{R) or / = f{R°'''Rab), there are, in 
general, two equations for the anisotropic stress in those models with / — f {R"'^'^'^ Rabcd) ■ Consequently, 
in addition to the isotropic and fully anisotropic metrics one must also consider the scenario of locally 
rotationally symmetric (LRS) solutions which without loss of generality are described by P2 = Pi = 
{H-p^)/2,H^3pi. 

For these LRS Kasner-like metrics, in the context of power-law Lagrangians, / = aZ", the field 
equations may be reduced to the two independent equations; 

= 8 - 8pi + 9pi^ - 4H - 6piH + H'' + 8n{H + pi - 2) (29) 

= ■iH^{n-l)-H'^{4-2?,n-8n) + 2H{l-2n){20-2^n + 'i2n) (30) 

+ (2n - 1)(47 - 148n + 128n^) - pi(n - 1)((5 - 8n)(13 - 16n) - i/(51 - 60n)) . 

For all n, these equations are solved by // = pi = 1, giving P2 = Pa ~ 0; this is the subset of the 
anisotropic solutions found in the previous section with H = J = 1 which corresponds to the Milne 
model. For n = 1 there are no other LRS solutions of this form. For n 7^ 1, the second of these equations 
is linear in pi; substituting the solution one obtains for pi into the first equation gives a sixth order 
polynomial equation for H in terms of n. However, factoring out the root H = 1 leaves the fifth order 
polynomial: 

= //''(n- l)^ + -H"''(l-n)(15-41n + 8n^)+-H"^(38-216n + 283n^ +56n^ -80n*) 
+H^{-M6 + 2156n - 5421n^ + 7032n^ - 4560n'' + 896n^) 
+H{1 - 2n)^(653 - 2594n + 3848n^ - 2688n^ + 1024n*) 
-(1 - 2n)^(361 - 2104n + 4768n^ - 4992n^ + 2048n*) . 

Consequently, in addition to the Milne solution, there must always be at least one locally rotationally 
symmetric Kasner-like solution with real values for pi and H for any real choice of n, other than n — 1. 
For example, when n — 1/2, there is a solution given hy pi = 1,H = 5, giving P2 = P3 = 2, whilst for 
n = 7/8 there are three real solutions, one of which is given by pi — —1/2, H = 1/2, so that p2 = P3 = 1/2. 
Thus it is dependent upon the choice of n as to whether the signs of the Kasner exponents must all be 
positive or not. These solutions are new and were not found in [25) . 

3.3.3 Isotropic power-law solutions 

In order to complete the discussion of this section, we consider the isotropic metrics, with pi = P2 = 
P3 = H/3. Except for Minkowski spacetime, which is a solution in any theory for which /(O) — 0, the 
relevant equation for these metrics in vacuum in the context of a power-law Lagrangian, / — aZ", is 

= -1 + 6n - 8n^ -f ^(1 - 2n + 4n^)H + |(n - 1)H^ . (31) 
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For n — 1, there is only one solution, with H = 3/2. For ti 7^ 1, there are two solutions, given by 

-l + 2n- 4^2 ± V-l + 10n-16n2 + 16n4) 



H = H± 



6(n~- 1) 



which are real if n < m = - (^5 + Vs + 2^5) /4 ^ -1.24 or n > n2 = (-V5 + Vs + 2^5) /4 « 0.12. 

Finally, we note that ii H ^ 3, equation (|3Hl is quadratic in n and so two values of n correspond to 
the same value of H and so Lagrangians which are a sum of two powers related by this equation also 
contain an isotropic power-law solution. Now H = 3 for n — i and thus, for / — OnZ" + OmZ"^ , with 



7^ 0, (4n— l)(n— 1) / and m = m^ 



with H = H^. 



5-16n + 20n^±3V-l + 10n-16u2+16n'' ,i ,i • i ,• 

„,. " ,., — TT , then there is a solution 

8(4n— l)(n — 1) ' 



3.3.4 Summary 

We have seen that anisotropic Kasner-like solutions in higher-order theories of gravity where the La- 
grangian depends only on the Kretschmann scalar, Z = RabcdR"'"^''^ , are possible only if the Lagrangian 
is of the form f{Z) = Z". Whilst it was claimed in an earlier work by Clifton and Barrow [25) that 
there are no exact anisotropic solutions of the form ([2)| in these theories, this is not quite correct. Firstly, 
for n = 1, the Kasner solutions of general relativity, for which H = J — 1, remain valid. This is to be 
expected, since the theory with / = aZ is a special case of quadratic gravity and all vacuum solutions of 
general relativity must also be solutions of the pure quadratic theory [33]. However, for all other values 
of n, fully anisotropic solutions of this form with the three Kasner exponents all real and distinct are 
not possible. Nevertheless, for choices of n other than 1, we have seen that solutions which have local 
rotational symmetry but exhibit an anisotropic third spatial coordinate are possible and the Kasner 
indices are subject to equations (|29|l and pop . 

In table [S] these results are summarised to show whether isotropic power-law, locally rotationally 
symmetric and fully anisotropic Kasner-like solutions are possible for various choices of the Lagrangian 
f{Z). Use is made of the definitions given in the previous section; ni = — ( \/5 + v 3 + 2\/5 ) /4 ~ —1.24 



and 712 = (-V5+ VS + 2%/5 ) /4 ^ 0.12 



fiZ) 




Minkowski 


Power-law 


LRS Kasner 


Anisotropic Kasner 


z 




/ 


/ 


X 


/ 


Z + A 




X 


X 


X 


X 


Z",n < ni 




X 


/ 


/ 


X 


Z",ni <n<0 




X 


X 


/ 


X 


Z",0 <n<n2 




/ 


X 


/ 


X 


Z",n>n2,nj^ 1 




/ 


/ 


/ 


X 


aZ" + PZ"'^ , n > 0, 


n/il 


/ 


/ 


X 


X 



Table 8: The existence of power-law and Kasner-like vacuum solutions in various models of f {Rabcd.R"''"^'^) 
theories of gravity. 



3.4 Quadratic gravity; f = KR + aW + PRabR"^ + iRabcdR""^'^ + A 

The theory of quadratic gravity, in which the Einstein-Hilbert Lagrangian of general relativity is supple- 
mented by quadratic Riemann, Ricci and scalar curvature corrections, is a particularly interesting special 
case to consider. The Lagrangian is given by 



f^>iR + aR^ + pRabR"'' + 'jRatcdR'''""' + A . 



(32) 



It was shown by Starobinsky [4] that addition of quadratic curvature corrections to the Einstein-Hilbert 
Lagrangian leads to the emergence of inflation. Furthermore, in contrast to general relativity, fourth- 
order gravity is renormalisable [ll] and thus it is often motivated as a first-order quantum correction to 
Einstein's theory. A review of the history of the study of these models of gravity may be found in [3]. 
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Without loss of generality, one may immediately set 7 = 0, since the Gauss-Bonnet term, defined by 
G = R^ — 4RabR"''' + RabcdR"^^'^'' , is a total divergence in four dimensions, so its variational derivative 
with respect to the metric does not contribute to the field equations. If the Ricci term is not present 
in the Lagrangian, i.e. /3 = 0, then the theory is a special of the f{R) theories studied earlier and, 
moreover, it is conformally equivalent to that of a minimally coupled scalar field in general relativity. 
Thus, one expects that the presence of the Ricci term in the Lagrangian might permit much more diverse 
anisotropic behaviour, as was found in [15]. The quadratic theory is scale- invariant iff kA — and 
conformally invariant iff k = A = and 3a + /3 — 0. 

In order to obtain the correct Newtonian limit in the slow- motion weak-field limit, it is necessary 
that the fourth-order terms contributed by the quadratic parts of the Lagrangian are exponentially 
vanishing, rather than oscillatory. To ensure that this is the case, the parameters must satisfy /3/k < 
and (3a 4- I3)/k > 0, with k / 0. 

3.4.1 Field equations 

For the line element given by ([2]| , the vacuum field equations for the theory of quadratic gravity defined 
by the quadratic Lagrangian in equation (|32|) are P^ = 0, and the relevant independent quantities are 
given by 

P: = ^ {{H^ ~ 2H + J)mH - 3){3a + 13) ~ Kt^) + AAt") (33) 



Po = ^{{3J-H'){3{J + H')a + {3J^3 + 2H)l3)+4H'{2H-3){3a + l3) 

+3{J - H^)Kt^ + QAf^) (34) 

P;-P; = ^(p^-p,.){{2{H-3){{H^-2H + J)a+[J~l)p)^{H-l)Kt') . (35) 

As before, ^, ;/ in equation ((35} for the anisotropic stress are not indices to be summed, but are used there 
as labels, taking values from 1 to 3. Otherwise summation convention is used as normal and overdots 
represent derivatives with respect to the coordinate time, t. 

One can immediately see that the cosmological constant. A, must be zero in order for any solutions 
of this sort to exist. 

3.4.2 K^O 

If the Einstein-Hilbert term is present in the Lagrangian, then the only possible isotropic power-law 
solution is Minkowski space, which requires A = 0. There is one family of anistropic solutions, that of 
general relativity, given by _ff = J = 1, provided A = 0. In light of the results of the previous sections, 
this is not unexpected. 

3.4.3 K = 

In pure quadratic theories, with k = 0, a non-zero energy-momentum tensor gives rise to a strong 
gravitational field and consequently spacetime is not asymptotically fiat [31], but we include them for 
completeness of this discussion. 

Any isotropic solutions of the field equations require A = and in general there are two such solutions; 
Minkowski space and the radiation-like solution with H — 3/2. For the special case of Weyl gravity, 
3a + P — 0, all power-law isotropic solutions are possible, since in the isotropic case, the contributions 
to the field equations from the _R^ and the Ricci terms in the Lagrangian are the same up to a constant 
multiple. 

There are several possible families of anisotropic solutions; as we have pointed out they all require 
the cosmological constant. A, to be zero. These solutions and the constraints that must be satisfied are 
summarised in table [9] 

Recall that the first class of solutions are vacuum solutions in general relativity and so have vanishing 
Ricci tensor. As a consequence, they must also be vacuum solutions of the quadratic theory for all values 
of the parameters a, /3 and k. Indeed, we have already seen that they are solutions for the special cases of 
quadratic Lagrangians within the more general models we have previously investigated. The constraint 
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Class Solution constraints 



Validity 



I 

II 
III 
IV 



H ^ J=l 

H = 3, 

J = st?(-3q + P- 2V-2/3(3q + /?)) 
H^3, 



This is a vacuum solution for all values of a, /3 
and K. 

These are solutions for the Weyl theory of 
gravity only, that is ii 3a + P — — k. 

These are valid with real exponents if «: = 
and -13/3 < a < -/3. 

These are valid with real exponents if «: = 
and -/3 < a < -/3/3. 



J = ST?(-3a + P + 2v/-2/3(3a + /?)) 

Table 9; Kasner-like vacuum solutions in pure quadratic gravity. 



on the solutions of class II implies that ps = pi + P2 — 1 ± 2-\/(pi — 1)(P2 — 1), and so all three Kasner 
exponents may be positive and greater than unity. Classes I-III were previously found by Deruelle [23] . 

3.5 Gauss-Bonnet theories; / = f{G),f = R + f{G) 

In four dimensions, the Gauss-Bonnet term is a topological invariant and its variation does not contribute 
to the field equations, though it may give rise to interesting cosmological effects in higher dimensions [35] . 
These terms arise in the low energy effective actions of string theory, however modified Gauss-Bonnet 
theories have also been proposed as a form of gravitational dark energy capable of successfully describing 
cosmology at late times [SS]. Here, we shall consider the situation where the Lagrangian is a general 
function of the Gauss-Bonnet invariant, i.e. / = f{G) and also the case f — R + f{G), where the 
Gauss-Bonnet term is defined by G = i?^ — ARabR"^ + RabcdR"''""^ , and takes the form 



G = ^{H - 3)(2pi' - 2Hpi +H^~J) 



(36) 



for the metric given by ^. 

The contributions P^ to the field equations due to a function /(G) in the Lagrangian are given by 



1 , . P1P2P3 



[{H - 3) fa - 3tfc 



-* 5. 



pM _ p^ _ 



-2/ + |(8pip2P3(^ - 3)fG + (J - H^){{H - 2)tfG + t^fa) , 



i3 



{p^ - p.) {H-p^~ p.) {H - 2) fa + tfa 



(37) 
(38) 
(39) 



where once again, fi, v in equation (]39|l are not indices to be summed, but labels taking values from 1 
to 3. Otherwise summation convention is used as normal. Here, /g is used to denote ^ and overdots 
represent derivatives with respect to the time coordinate t. 

The anisotropic stress will vanish independently of / both for isotropic solutions, ie pi = P2 = Ps, and 
solutions in which two of the Kasner exponents are zero, and also in theories in which the Lagrangian is 
a power of the Gauss-Bonnet term, provided that H is suitably chosen in this case. 

3.5.1 / = /(G) 

Except for the trivial case of /(G) = G, the only possible real anisotropic solutions have G = 0. Fur- 
thermore, for any Kasner solutions to be possible, it is required that /(O) = 0. Provided /g(0) does 
not diverge, then any solution of G = defines a two-parameter family of exact Kasner-like solutions. 
If, on the other hand, /g(0) is divergent, then for general such /(G), there is a one-parameter family 
of solutions where the three Kasner indices are given by pi = P2 = 0, with the third index free, plus 
permuations. In the special case /(G) = aGlogG, there are additional families of solutions given by 
Pi — 0, with p2 and pz free (and permutations), or J = 2pi^ — 2_H'pi -I- H^ . 
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3.5.2 f = R + f{G) 

The only real solutions with G — and Rab = are Minkowski space and the Milne model, corre- 
sponding to pi = l,p2 = P3 = (plus permutations), which are vacuum solutions whenever /(O) = 0. 
Consequently, only these solutions can separately satisfy both the vacuum Einstein equations. Gab ~ 0, 
and the vacuum equations due to the purely Gauss-Bonnet terms, Pab ~ 0. Other solutions must have 
G 7^ and the Einstein-Hilbert terms must balance the Gauss-Bonnet terms, that is to say they must 
be of the same order in time. Therefore, the possible real Kasner-like vacuum solutions in these theories 
may be summarised as in table [10] 



Class 



Solution constraints 



Validity 



This is a vacuum solution if / = aVG + 

^ log G. 

This is a vacuum solution if / = aG log G— 

3V2G- 

This is a vacuum solution if / = a\fG. 

These are vacuum solutions if / = a\fG. 

These are vacuum solutions if / = a\fG. 
These are vacuum solutions if / = aVG. 



I 

II 

III 
IV 

V 

VI 



n — P2 — P3 



Pi 



- P2=PZ 



Pi^P2= P3 = 63^(3 + a^± \/12q2 + 9) 



Pi = 1,P2 =P3 = 1 +«^ 

(and permutations) 



H = 1,J = l + 4pia^ +Apia^l - pi + a^ 



H = 1,J = l + 4pia^ -Apia^l - pi + a^ 



Table 10: Povifer-law and Kasner-like vacuum solutions in a modified Gauss-Bonnet gravity. 

We can see that although there are several different and interesting types of Kasner-like solutions in 
these theories, they are only valid for a small set of Lagrangians. This is because only the Minkowski 
and Milne universes solve both the vacuum Einstein equations, Gab = 0, and the vacuum field equations 
due to the purely Gauss-Bonnet terms in the action. Thus, terms in the field equations due to the 
Einstein-Hilbert term must be balanced by those due to /(G) and so these terms must be of the same 
order in time. 

3.6 Weyl theories of gravity 

We now consider the situation where the Lagrangian is a general function of the Weyl invariant, i.e. 
/ = f{W), where the Weyl term is defined by W = \R^ - 2RabR'''' + RabcdR"'""^ ■ We may write 



W" = -^ ((Pl -P3)(P1 -P2)(pi - 1) - (P2-P3)(pi -P2)(P2- 1) + (pi -P3)(P2 -p3)(p3 - 1) ) 

If we assume, without loss of generality, pi > P2 > P3, then we see that the first and third terms are 
positive, but the second is negative. However, the middle term must be no greater in absolute magnitude 
than the third term if p2 < 1 and similarly it must be no greater in magnitude than the first term if 
P2 > 1. Thus the Weyl tensor is non-negative and is zero only if each term vanishes separately. This 
requires that either the metric is locally rotationally symmetric with wlog p2 = P3 and pi = 1, or it is 
isotropic, pi = P2 = P3. The vacuum field equations are P^ = 0, with 

PI = ~2f + 2Wfw, (40) 

PS = -^f + ^{lipiP2Pi + H^-2H''-iJ + AHJ){{H-i)fw+tfw)). (41) 

PH -P: = ^ {P^-P.) (l2t{H^p^-p,)i(H-2)fw+tfw) + 4{H ~ 3)(2(J- 1) ~ {H - lf)fw 



+2{3{H - 1) - 5{H - if + 4(J - l))tfw - 6{H - l)t\fw 



(42) 



where once again, fj,, v in equation (|39p are not indices to be summed, but labels taking values from 1 to 
3. Otherwise summation convention is used as normal, fw is used to denote f'(W) = -^ and overdots 
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represent derivatives with respect to the time coordinate t. The equations (|42|) for the anisotropic stress 
may be combined to obtain 



= ^ (pi - P2) (pi - Pb) (P2 - P3) ({H - 2)/w- + tfw) 



(43) 



For fully anisotropic solutions, the Weyl tensor is non-zero, and as before one can treat the field 
equations as differential equations for /, solving them to find the required form of Lagrangian for the 
higher-order theory to possess such solutions. It is found that the only possible choice is that of Weyl 
gravity, f{W) = aW, and that solutions require pa = pi -I- P2 — 1 ± 2^{pi — l)(p2 — 1); this is a special 
case of the quadratic Lagrangians considered before in section \3A\ corresponding to 3a + p = 0. 

For solutions of the vacuum field equations where the Weyl tensor vanishes, it is necessary that 
/(O) = 0. In fact, if /(O) = and f'{W) does not diverge at W = 0, then all solutions oi W ^ are 
solutions of the vacuum field equations. However, if /(O) = and /'(O) diverges then we must consider 
both the isotropic and locally rotationally symmetric types of solution in turn as the existence of such 
solutions will depend upon the manner of this divergence. 

These results are summarised in table [TT] 



Class 



Solution constraints 



Validity 



I 



Ila 



lib 



III 



IV 



Pi = P2 = Pi 



Pi = P2 = P3 = 2n - 1 



Pi = P2 = P3 = (4n - l)/3 



Pi = 1,P2 =P3 

(and permutations) 

Pi = 1,P2 =P3 = 2n - 1 
(and permutations) 



P3 = pi -I- P2 — 1 ± 2 a/ (pi — l)(p2 — 1) 



These are vacuum solutions for all values of pi if 
/(O) = and either /'(O) converges or diverges 
due to terms in f{W) of the form W" with 1/2 < 
n < 1. 

This subset of the class I solutions remains a vac- 
uum solution if /(O) = and the divergence in 
/'(O) is due to terms in f{W) of the form W" 
with < n < 1/2. 

This subset of the class I solutions remains a vac- 
uum solution if /(O) = and the divergence in 
/'(O) is due to terms in f{W) of the form W" 
with < n < 1/2. 

These are vacuum solutions if /(O) = and /'(O) 
either converges or diverges due to terms in f{W) 
of the form W" with 1/2 < n < 1. 

This subset of the class III solutions remains a 
vacuum solution if /(O) — and the divergence 
in /'(O) is due to terms in f{W) of the form VK" 
with < n < 1/2. 

This is a vacuum solution in Weyl gravity, f{W) = 
aW. 



Table 11: Power- law and Kasner-like vacuum solutions in Weyl gravity. 

Note that if /(W) contains a term VW, i.e. n — 1/2, then the only solution is that of the Milne 
model, with metric given by 

ds = —dt +tdx + dy + dz , 

which we have seen is Minkowski space in a different coordinate system. Finally we point out that if the 
f'{W) diverges at W = due to a type of term other than a power law, there are no vacuum Kasner-like 
solutions in that theory. 



3.7 Homogeneous Lagrangians 

Thus far, this study has considered models in which the Lagrangian is a general function of one of the 



curvature scalars, X = R, Y 



K'^'Rab, or Z 



R'^ '^ Rated, or of a particular combination of these. 



17 



and <f) = 


^, and further define 




da 
86 ' 




(45) 


da 

d(t> ■ 




(46) 



Consequently, in each case there has been only one timescale in the problem, and, by substituting for 
the time variable using the appropriate curvature invariant, it has proved possible to derive all possible 
solutions and their existence conditions within these wide classes of models. 

However, it is useful to also consider more general Lagrangians depending on more than just one 
variable. We may also conisder Lagrangians which are homogeneous in X^ ,Y and Z, that is to say for 
all values of A, /(AX^, AF, \Z) = A"/(X^, Y, Z) for some n. X^ is chosen as the argument here, rather 
than X, in order that each term in the function will be of the same order in time. Thus, functions of 
this form will be relevant where the dominant terms at early or late times are more complicated than in 
those Lagrangians studied previously, such as a monomial in X, Y, Z. We may write 

f{X,Y,Z)=X^''a{e,<l>), (44) 

where a is a general differentiable function of 6 -. 

7 

It is important to note that a, /3 and 7 take constant values which are dependent upon the exact form 
of the Lagrangian / and also the values of 9 and </> for a particular solution. For example, in the case of 
the monomial given by 

f ^S^X^Y^Z", (47) 

we have 

2n = cT + 2n + 2v, (48) 

" - K^) (I) • <") 

l> = -(t> + -)|j", (50) 

7 = My"- (51) 

Inserting the general form (|44p for the homogeneous Lagrangian into the vacuum field equations for 
the metric ((2|, one obtains the equation 

= —^^ — (pi - P2) (pi - pa) (P2 - pa) {n-l){l + H- 4n) (/3X' + 7^) . (52) 

We will consider in turn the solutions of this equation in the remaining field equations. 

3.7.1 Isotropic solutions 

In the isotropic case, pi = P2 ~ Pa, there is only one independent field equation, given by 

= pi(pi(2pi-l))'""' (a(2p?-2pi + l)'(l + 2pi(n-l)-6n + 8n^) 
-(6/? + 7)pi(l - 2pi)2(l + 3pi - 4n)) . 

Recall that a, j3 and 7 are constants which will depend upon pi according to the exact form of the 
Lagrangian /. Thus, given a choice of /, this equation provides a necessary and sufficient condition for 
there to be a isotropic power-law solution expanding as t^'- . 
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3.7.2 Locally rotationally symmetric solutions 

In the case of locally rotationally symmetric (LRS) solutions, there are two independent field equations. 
Without loss of generality, we set P2 = Pa 7^ Pi ■ 

LRS solutions with X = are a special case of the fully anisotropic ones of that kind and will satisfy 
the same existence conditions. Similarly, if pi + 2p2 = 4n — 1, then there are LRS solutions provided 
aX^" = can be solved simultaneously. This requires either the additional constraint a = to be 
satisfied, or if 1/4 < n < 5/8 then there are real solutions of this form with X — 0. Both these types of 
LRS solution satisfy the same constraints as in the fully anisotropic case, and so these will be discussed 
at greater length in the next section. 

If we define X = t^X, Y = t^Y and Z = t'^'Z for convenience, then, for X 7^ 0, these may be reduced 
to 



= InaZ' -2X{2-3pi-2p2 + {pi~P2f+in{Pi+P2-^))il3X''+-fY)+l3X''Z 

+-fXZ{pi' + 2pi - 1 + 2(pi + 2p2 - l)(n - 1)) j (1 +pi + 2p2 - 4n) , 
= aZ" (^4n(pi^ + 2p2'' + (4n-3)(pi+2p2)) -l) 

+ 16Xp2'(l + pi + 2p2 - 4n) (pX {pi' - 2p^p2 + P1P2' - ^Pi - 2piP2 - 2p2'') + 7((2 - 3p2)p2^ 
+pi' + 2pi(p2 - 1)P2 - 2pi*(l + P2) + Pl^P2(6 - 7p2 - 4p2^) + p/(l - 4p2 + 3p2')) 



In this way, the field equations give relationships between the theory-dependent constants a, /3 and 7 
which must be satisfied if the theory is to contain any other LRS solutions. 

3.7.3 Anisotropic solutions 

We now consider those solutions of equation (|52|l which are fully anisotropic with pa all distinct. We 
have seen that for metrics described by the line element ^, Z = Q only for Minkowski space and the 
Milne model. Thus for fully anisotropic solutions, Z will be non-zero. The only real Kasner-like solutions 
with X = Y = Q are those of general relativity, which satisfy H = J = 1. Since the leading order terms 
in the field equations as X, y — >■ behave like a{2n — l)X^"~^ , then if a diverges at X = 0, faster than 
■^i-2n^ there will be no solutions. However, if a is finite at F = 0, one either requires that n > 1/2 or 
that a — and n > so that the next leading order terms do not diverge. For solutions with X = but 
Y y^ 0, the leading order terms in the field equations as X — >■ behave like a{H — in + 1)X^"~^, so if a 
remains finite at X = 0, one either requires that n > 1/2 or that either H = 4n~loTa = and n > 
so that the next leading order terms do not diverge. 

If there are configurations of the Kasner exponents such that a — 13 — y = 0, then these will be 
solutions of the theory. Recall that a, j3 and 7 depend on the explicit form of the Lagrangian. 

Substituting // = 4n — 1 into the field equations leads to the necessary and sufficient condition 
aX^" _ Q Pqj. j^ _ q -^vith _H" = 4n — 1, this implies that J = 16n(l — n) — 3 and requires that 
1/4 < n < 5/8 for these solutions to be real. Again, since a depends on the explicit form of the 
Lagrangian and the values of the Kasner exponents and so a = defines an extra constraint on the 
solution. 

In addition to those solutions which are a special case of the ones with H = An — 1 and a = Q, there 
is another class of solutions if n = 1. This class of solutions exists if 

a = __L(x^_4y- + ^) and (53) 

P = -^(4y-Z) . (54) 

Again, a, /3 and 7 will depend on the values of the Kasner exponents through the curvature scalars X, Y 
and Z and so these equations will typically give two further constraints on the possible set of solutions. 
However, if for example, a — a(X^/(4y — Z)), then the second of these equations will be trivial and if 
a = A(X^ — 4y -|- Z)/X^ for some constant A, then both of these will be trivial. 
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Finally, there are several additional solutions if j3X'^ + ^Y = 0. This equation is satisfied trivially 
if the Lagrangian / is independent of Z. Otherwise, it provides a constraint on the solutions. The 
remaining field equations then show that if n = 1 there are solutions subject Xo H = 1 and 

«=|^(1-^)- (55) 



If n = 1/2, there are solutions subject to 



X^ 



^ = '2P^{1-H){H^-J). (56) 



Finally, for general n there are solutions of this sort if 

a = 2P ^ {H-l){H-H^ + n{H^-2H + J)) and (57) 

J = —{l-2H + H'^-6n + 8Hn + ?,n - 9,Hn) . (58) 

3.7.4 An explicit example 

As an explicit example, we shall consider fully anisotropic Kasner-like solutions in the model of Briining, 
Coule and Xu 37] , which is described by the Lagrangian 

/ = ^ + Aj+r|, (59) 

where A and r are constants, and was studied there in the context of FRW solutions. Note that the 
model reduces to general relativity in the limit A = r = 0, so in the discussion that follows we will assume 
that A and r are not both zero. In terms of our notation, 



n = 


2 ' 


a - 




P = 


-- -^{\Y + rZ) 


7 = 


^ ^^^ 



and 



Thus we can explicitly observe the dependence of a, P and 7 upon the values of the Kasner parameters pa 
and in particular for this model we note that a diverges for the general relativistic solution, H = J = 1. 
For anisotropic solutions in this model, with the Kasner exponents pa all distinct, equation (|52l) 
reduces to 

O-rf-Ji-^) , (60) 



^H2-2H + J^ 

and so the possible solutions depend on the value of r. For the special case when r is zero, the field 
equations reduce to the single equation 



2H^ - 3H'^ + J + 2HJ - H^J - J= 



{H^-2H + jy 



A = 1 . (61) 



This equation allows a two-parameter family of solutions to be found with real- valued Kasner exponents 
for all values of A except those in the range — 1 < A < 0. 

For the general case with r non-zero, solutions require H = 1 and 

C(J - if + 8pi(J - 1) - 16pi'(pi - 1) = , (62) 

where we have defined ^ = (1 -I- 3 A -I- t)/t. These solutions therefore are described by the single free 
paramter pi. Seeking solutions for which the Kasner exponents po, take real values, it is necessary to 
choose pi such that 

l+C(Pi-l)>0 (63) 

20 



and either 



ipi + 2 ^i - 1 j pi - 1 < or (64) 

?(C(l + 3pi)'-16pi) < 0. (65) 



It is possible to find pi such that the inequalities (|62p and (|64p may be satisfied simultaneously if C < 5/4, 
whilst the inequalities (|62|) and (|65p may be satisfied simultaneously provided ^ < 4/3. Thus, this bound 
provides the restriction on the allowable choices of A and r such that the general model (|59|) will contain 
anisotropic Kasner-like solutions of the form p]). 

If we compare these findings with those in the previous section, when the general case of a homoge- 
neous Lagrangian was discussed, we see that these correspond to the solutions of equation (|56p . When 
r = 0, the equation /3X^ + 7F = is trivial, but for non-zero r, it provides an additional constraint 
and reduces the number of free parameters in the solution from two to one. There are no anisotropic 
solutions corresponding to X = 0, since (except in the limiting case of general relativity) aX diverges 
there. 

3.7.5 Remarks 

In this section, we have studied those Lagrangians which are homogeneous in time for the metrics Q, 
and which may be written in the form (|44p . We have found the existence conditions for all possible 
Kasner-like solutions in this class of theories. Although it is impossible to find the exact forms of these 
solutions without explicitly defining the Lagrangian, and the dependence of the functionals a, /3 and 7 
on the Kasner parameters may be complicated, we have developed a general framework which can be 
used to obtain all possible solutions once the gravitational theory has been defined. This framework was 
then explicitly demonstrated for the model of Briining, Coule and Xu [3Z| • 

Since the terms in the field equations must vanish at each order in time, then unless there exists 
a particular configuration of Kasner exponents such that the Lagrangian f[X, Y, Z) is a constant, it is 
reasonable to expect that the solutions found in this section will be the only Kasner-like solutions with 
real Kasner indices that are possible in higher-order metric theories of gravity derived from Lagrangians 
with the general form f{R, RabR"'', Rated R"'""' ) ■ 

4 Anisotropically Inflating Solutions 

In the preceeding sections, we have studied higher-order theories of gravity in which the Lagrangian is 
dependent upon the scalar curvature, X = R, the Ricci invariant, Y = RabR'^^ , and the Kretschmann 
scalar, Z = RabcdR°'^'^^\ to find the conditions required on the form of the Lagrangian for a particular 
theory to contain Kasner-like vacuum solutions, where the line element is described by ^. Within the 
context of this class of higher-order gravitational theories, it is interesting to also study the possibility 
of other simple anisotropic exact solutions. A natural yet simple extension is to consider a Bianchi type 
I solution which describes exponential but anisotropic expansion, where the metric is of an anisotropic 
deSitter-like form; 

ds^ = -dt^ + e^^^'dx" + e^'^^'dy^ + e^^''dz\ (66) 

and the exponents pa are real. 

In this section, we shall investigate whether such metrics solve the vacuum field equations of these 
theories and we include also the possibility of a cosmological constant. Such solutions were found to 
be possible in theories with quadratic Lagrangians, provided that the Ricci term is present [15] and 
their existence demonstrates that the cosmic no-hair theorem of general relativity, which states that the 
presence of a positive cosmological constant drives the solution towards the deSitter one at late times, 
cannot be extended to higher-order theories in general. 

The three curvature scalars X, Y, Z are constants, and take the values 

X = H^ + J , 
Y = {H^ + J)J, 
Z = 2(J^ + M) , 
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where the useful definitions H = pi + P2 + Ps, J = Pi + P2 + P'i and M = pi + p2* + P3 have 
been made. Note that all three curvature scalars X, Y and Z are constant and also positive unless 
pi = P2 = Pi = 0, corresponding to Minkowski space. Thus, unlike the situation for the Kasner solutions 
discussed previously, it is not necessary for example that R = Q for f{R) to be zero. Furthermore, since 
the curvature scalars are constant, the field equations simplify substantially [3H], and for the isotropic 
case they reduce to the single equation 

\f=^f.^lH^fr + ^^H^fz, (67) 

which must be satisfied by / if the de Sitter universe is to be a solution in a particular higher-order 
theory of gravity of the form /(X, Y, Z). 

For all vacuum anisotropic solutions of the form (|66p . they become 

/ = 8Hpip2P3fz, (68) 

fx +2.7 fy = 2iH^-3J)fz, (69) 

where, as before, subscripts are used to denote differentiation of / with respect to that curvature scalar, 
so that, for example, fx = w- Any cosmological constant is to be included in the function /. 

If the Lagrangian is a function of only one of the three curvature invariants, / — /(^) say, where 
^ £ {X, Y, Z}, then the field equations for anisotropic solutions simplify further to 

/ = 0, (70) 

/« = 0. (71) 

Thus the theory will contain a two-parameter family of exact anisotropic deSitter-like solutions if the 
Lagrangian, /, has a positive double root, ^ = ^o > 0, the simplest example of which being a quadratic 
Lagrangian, / = a{£, — ^o)^. It should be noted that for the case where ^ = R, the scalar curvature, 
choosing the sign of the Lagrangian to give the correct Newtonian limit would give a negative cosmological 
constant due to the required positivity of Ro- 

In this way, some degree of fine-tuning is required in the Lagrangian for such solutions to exist. 
However, more general types of model do not require such fine-tuning for solutions of this sort to exist. 
In particular, if the Lagrangian is of the form / — f{Y/X^), then equation (I69|) is identically satisfied 
and so it is only required that there exist positive roots of / = 0. 

5 Conclusions 

We have studied some aniosotropic cosmological Bianchi type I solutions to a wide class of higher-order 
theories of gravity derived from the three curvature invariants R, RabR'^'' and Rated R"'"^'' ■ Although 
general relativity is well-supported by solar system tests, in the high curvature limit, such as on approach 
to an initial cosmological singularity, we expect quantum efi'ects to become important and to cause 
deviations from the standard behaviour in general relativity. At high curvatures, anisotropics diverge 
faster than isotropics and will tend to dominate the cosmological behaviour at early times. Furthermore, 
it has previously been shown [T5] that anisotropies in these higher-order theories may display significantly 
different behaviour to that found in general relativity. Thus, in this work we have investigated the role 
of anisotropy on approach to the initial singularity. 

In particular, we have found all Kasner-like solutions given by the Bianchi type I line element ([SJ 
for several wide classes of higher-order theories of gravity, and also all of the similar Bianchi I solutions 
which are described by the line element (|66p . Previously, Kasner-like solutions were known only for 
quadratic gravity (25] and higher-order Lagrangians which were powers of one of the curvature invariants 
RjRabR"''' and RabcdR"'"^''^ ■ We have extended this to much more general Lagrangians including those of 
the form ,f (-R), /{RatR"^'') and f {RabcdR'^^'^'^) , and in the course of this investigation, we have also found 
additional solutions to the previously-studied power- law Lagrangians which were not found in 24ii25]. 

We have further widened this study to include also those Lagrangians which are homogeneous in time 
for the metrics ((5}, which may be written in the form (|44p . using the model of Briining, Coule and Xu [3Z] 
as a particular example. Since the terms in the field equations must vanish at each order in time, then 
unless there exists a particular configuration of Kasner exponents such that the Lagrangian f{X, Y, Z) is 
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a constant, we expect that the solutions found in section 13771 will be the only Kasner-like solutions with 
real Kasner indices that are possible in higher-order metric theories of gravity derived from Lagrangians 
with the general form f{R, RabR"'', RabcdR"'""^) ■ 

Although they are geometrically special, the Kasner-like solutions given by the Bianchi type I line 
element © provide us with a very useful insight into the dynamics of anisotropics, and also they give a 
good description of the evolution of more general anisotropic cosmological models over finite time inter- 
vals. This is of particular interest when considering the behaviour on approach to the initial singularity 
exhibited by the Bianchi type VIII and type IX ( "Mixmaster" ) cosmologies, which can be approximated 
by a sequence of different Kasner epochs. We have considered the properties of the solutions found in 
relation to the behaviour of these more general anisotropic cosmologies and found that in general it is 
model-dependent as to whether the universe will experience an infinite sequence of oscillations between 
Kasner regimes as the singularity is approached. A more detailed analysis is required to understand the 
extent of the validity of these vacuum solutions in the presence of a non-comoving perfect fluid. 

The conditions for the existence of de Sitter solutions in higher-order theories of gravity were already 
known |38) . and some examples of anistropically inflating solutions had previously been discovered for 
particular theories [15]. However, we have extended this work and explicitly discovered the existence 
conditions for anistropically inflating solutions of the form (|66p in all higher-order metric theories of 
gravity derived from Lagrangians with the general form f{R, RatR"^^, RabcdR'^'"^'')- For the simpler models 
in which the Lagrangian depends on just one of the curvature invariants R, RabR"''', and RabcdR"''"^'^ , we 
found that some degree of fine-tuning is required in the Lagrangian for such solutions to exist, however 
this is not required in more general theories. These solutions explicitly demonstrate that the cosmic 
"no-hair" theorem of general relativity does not hold in general in higher-order theories. 
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